General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



Laboratory for Information and Decision Systems 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 


(HASA-CB-158474) SUSOLiTIOM BTALUATI09 OF 
C0HBI9ED 4D BEAT AWD I If BOB BE TRAFFIC 
SIT0ATI0B DISPLAYS ABD PROCEDURES APPLIED TO 
TBBHIHAL AEBIAL HAW02YEBS Semi-annual 
Progress Report, 1 Sep. 1978 (Bassachvsetts G3/04 


■ 79-21033 


UFnclas 

1*735 


Semi-Annual Progress Report 


Simulation Evaluation of Cochined 4D RNAV and Airborne 
Traffic Situation Displays and Procedures Applied 
to Terminal Area Manuevers 


Principal Investigators: Michael Athans 

Mark E. Connelly 


Period: September 1, 1978 to March 1, 1979 







(NASA Ames Grant NSU-2180) 


April 1979 




Introduction 


initial experiments show that the aircraft arrives at the 3D waypoint within a 


few seconds of the anticipated time. The presence of headwinds or tailwinds 


of 


The 






been written by the project based on 

' 

the descent algorithm analyzed in the last semi-annual progress report: 

1. Complete Algorithm for the Ti59 Hand Calculator (predicts hori- 
zontal distance and elapsed time for an idle thrust descent from 
36,000’ to 10,000’ given the desired Mach during the first phase of 
the descent and the desired equivalent airspeed during the second 


phase). Running Time is 5 minutes, 13 seconds. 
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is completed. As a consequence, the computation must be repeated (iterated) 
one or more times to obtain more accurate approximations for the missing values. 
In the M.l.T. descent algorithm, no more than two passes through the routine at 
each altitude level are employed. Additional iterations, it was found, have 
very little effect on the results. 





The same basic algorithm is implemented quite differently in FORTRAN, in 
the TI59, and in th* Adage assembly language ADEPT because of the unique con- 
straints imposed in each case. The FORTRAN program, for example, makes one 
complete pass through all the alcitude levels, storing the iterim results in 
tables for use in the second pass. This extravagant use of storage tables 
was impossible in the T159, which is memory limited. In the latter program, 
consequently, two iterations are executed at each altitude in sequence in the 

full version and one iteration per altitude in the short version. Minimizing 



the running time was the main consideration in organizing the real-time ADEPT 



program. Since it operates repetitively at a one solution per second rate 
during the descent, values generated in the prior frame are used to approximate 
the missing variables in the current frame, hence only one pass through the 
algorithm is required per frame. 

A flow chart for the FORTRAN version of the algorithm is given in Figure 1 
and a flow chart of the short TI59 program is presented in Figure 2. A listing 
of the FORTRAN program is given in Appendix A as well as a set of outputs for 
three descent profiles (standard, medium speed, slow speed) . A listing of the 
short TI59 program is given in Appendix B. The real-time Adage descent program 
is still being tested and refined, hence it will be documented at a later date. 

The results from the FORTRAN and TI59 program agree quite well. For example. 





on a standard .93 Kach - 320 knot EAS descent from 36,000' to 10,000* the FORTRAN 
program predicts an elapsed time of 591 seconds and the TI59 program predicts 
591.6 seconds. Horizontal distance according to the FORTRAN program will be 
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Figure 2 
15-7 New sHour fti gokithm 
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439,996 feet, whereas the TT.59 estimates 441, 165 feet 
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Two visiting students from the Centre d* Etudes et de Recherches de Toulouse 
M. Roger Dubois and M. Jean-Michel Loussant, will be working full time from 
April through June on the evaluation and refinement of the TI59 algorithm. 


707 Simulation Proqram Refinements 


The decision to include open-loop descent tests in the research program 

was based on the need for an interim 4D RNAV technique that could be implemented 

in the near future. Open-loop testing, however, places much more severe re- 

... 

quirements on the fidelity of the simulation model. Modeling errors that would 
be washed out in a closed-loop system, cause cumulative performance errors in 
an open- loop system that would completely mask the human factor effects we are 
attempting to measure. For this reason, a very significant effort has been 
made to correct errors and to improve the accuracy of the M.I.T. 707 simulation 
prooram. The original 707 program written by Captain Charles Corley (reference 
1) has been thoroughly checked and over a hundred major modifications made to 
enhance its fidelity. A separate report on the current aircraft model and 
simulation program is now in preparation, so we will not go into great detail 
here. A structural flow chart for the new 707 program is given in Fig. 3. 

The following list identifies the most important changes: 

1. The roll moment equation was modified to include the inter- 
action between inboard and outboard ailerons as a function of flap 
setting and the effects of spoiler blowdown. 

2. The airspeed dial now reads equivalent airspeed (EAS) in- 
i’, tead of indicated airspeed (IAS) since the profile descents utilize 
constant EAS, which is a true indicator of dynamic pressure. 


3. The alignment of the MACH dial and the EAS dial is now cor- 


rect at the current airspeed needle position. The MACH dial range 
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was reduced to .6 through .9. 

4. The rotating drum on the altimeter (100's of feet) was re- 
placed by an discretely updated digit. 

5. The atmospheric variables and the instrument calculations 

fa 

ate now updated every 1/15 seconds. 


6. All the analog control inputs have been rescaled as well as 
the corresponding aerodynamic coefficients - 

7. Errors in drag, the Euler angle routine SICGS, thrust, 
elevator limits, and several logic errors have been corrected. 

8. The old piecewise-lxnear traffic generator has been removed 
to make way for the new traffic routine which will be based on stored 
trajectories. 

9. Three new function switch modes have been created: 

(a) Fly with position frozen (useful in establishing the 

steady-state flight conditions required by the initial 


condition routine INITK) . 


(b) Return to monitor leaving the real-time program in a 
clean condition so that it can be restarted at the point 
where the interruption occured. 

(c) Execution of a single frame (used for debugging) . 

10. The steady-state elevator bias can now be set as an initial 
condition, thus avoiding a severe pitch transient when the program 
is started. 


11. A more accurate square root routine has been written and the 
duplication of square root software in FNST and ACSM eliminated. 

12. The trapezoidal integration of acceleration components has 


been rescaled to reduce the risk of overflow. 





M.S. Thasis, M.I.T. Department of Electrical Engineering, December 1974. 
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